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ABSTRACT 

We present the first X-ray monitoring observations of the X-ray bright FR I radio galaxy NGC 6251 
observed with RXTE for 1 year. The primary goal of this study is to shed light on the origin of the 
X-rays, by investigating the spectral variability with model-independent methods coupled with time- 
resolved and flux-selected spectroscopy. The main results can be summarized as follows: 1) Throughout 
the monitoring campaign, NGC 6251 was in relatively high-flux state with an average 2-10 keV absorbed 
flux of the order of 4.5x 10~^^ erg cm^^ s^^ and a corresponding intrinsic luminosity of 6x10^^^ erg s~^. 

2) The flux persistently changed with fluctuations of the order of ~2 on time scales of 20-30 days. 

3) When the hardness ratio is plotted against the average count rate, there is evidence for a spectral 
hardening as the source brightens; this finding is confirmed by a flux-selected spectral analysis. 4) The 
fractional variability appears to be more pronounced in the hard energy band (5-12 keV) than in the soft 
one (2.5-5 keV). 5) 2-month averaged and flux- limited energy spectra are adequately fitted by a power 
law. A Fe Ka line is never statistically required, although the presence of a strong iron line cannot be 
ruled out, due to the high upper limits on the line equivalent width. The inconsistency of the spectral 
variability behavior of NGC 6251 with the typical trend observed in Seyfert galaxies and the similarity 
with blazars lead support to a jet-dominated scenario during the RXTE monitoring campaign. However, 
a possible contribution from a disk-corona system cannot be ruled out. 

Subject headings: Galaxies: active - Galaxies: jets - Galaxies: nuclei - X-rays: galaxies 



1. INTRODUCTION 

Non-blazar radio-loud active galactic nuclei (AGNs) - 
objects with jets forming large viewing angles to the line of 
sight- are traditionally divided into two classes: Fanaroff- 
Riley II (FR II), with 178 MHz powers > 2 x lO^^ W/Hz 
and edge-darkened radio morphologies, and FR I galaxies 
with lower powers and more compact morphologies (Fa- 
naroff & Riley 1974). For the same host galaxy optical 
magnitude, FRIs produce about one order of magnitude 
less optical line emission than FRIIs (Baum et al. 1995) 
and have fainter or negligible UV continuum fluxes (Zirbel 
& Baum 1995, Ho 1999). 

While earlier models for the origin of the FRI/II di- 
chotomy focused mainly on accounting for their large-scale 
radio morphologies, more recently new ideas have emerged 
concerning the nature of the central engine in the two types 
of radio galaxies in an attempt to explain the nuclear prop- 
erties. One school of thought is that the nuclear X-ray 
properties of FR I and FR II are related to a different 
accretion rate onto the central supermassive black hole, 
with FRII being dominated by relatively large values of 
L/LEdd, while FRI would be accreting at sub-Eddington 
rates, L/LEdd ^ 10"'^ (e.g., Reynolds et al. 1996; Ghis- 
ellini & Celotti 2001). 

A fundamental step to gain insight into the nature of 
the central engine in radio-loud AGNs is to understand 
whether the X-ray radiation is produced by disk/corona 
systems as in Seyfert galaxies or by jets as in blazars. 



While the time-averaged spectroscopy (due to spectral de- 
generacy) and the pure temporal analysis (due to the fact 
that both radio-quiet and radio-loud show strong variabil- 
ity) cannot firmly discriminate between the two competing 
scenarios, their combination, i.e., time-resolved spectral 
analysis and energy-selected temporal analysis, offers in 
principle a better way to distinguish between accretion- 
dominated and jet-dominated systems. This conclusion is 
supported by the strikingly different spectral variability 
behavior shown by Seyfert-like objects (e.g., Papadakis et 
al. 2002; Markowitz & Edelson 2001) and by blazars (e.g., 
Zhang et al. 1999; Fossati et al. 2000; Ghozzi et al. 2006). 

Here, we concentrate on the X-ray nuclear properties of 
the nearby radio galaxy NGC 6251 (z=0.024), which is a 
giant elliptical galaxy hosting a supermassive black hole 
with mass Mbh ~ 4 - 8 x lO^M© (Ferrarese & Ford 1999). 
Based on its radio power at 178 MHz, NGC 6251 is clas- 
sified as an FR I (e.g., Owen & Laing 1989), whereas in 
the optical, it is classified as a type-2 AGNs (e.g., Shuder 
& Osterbrock 1981). 

Despite the intensive study of this source at all wave- 
lengths, the nature of the accretion process in NGC 6251 
is still a matter of debate. Based on the radio-to- X-ray 
spectral energy distribution. Ho (1999) suggested that an 
Advection-Dominated Accretion Flow (ADAF) is present 
in the nucleus of NGC 6251. On the other hand, Fer- 
rarese & Ford (1999) and Meha et al. (2002) favored a 
standard accretion disk. Finally, Mukherjee et al. (2002) 
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and Chiaberge et al. (2003) advocate a jet origin for the 
broad-band emission, based on the possible association of 
NGC 6251 with the EGRET source 3EG J1621+8203 and 
on its spectral energy distribution, respectively. 

In the X-ray band, NGC 6251 has been previously 
observed with various satellites. For example, ROSAT 
showed the presence of an unresolved nuclear source em- 
bedded in a diffuse thermal emission with temperature 
kT ~ 0.9 keV, associated to the galaxy's halo (Birkin- 
shaw & Worrall 1993). In addition, the existence of a 
correlation between soft X-rays and radio fluxes prompted 
Hardcastle & Worrall (2000) to hypothesize a jet origin 
for the soft X-rays. On the other hand, the prominent Fc 
Ka line (EW~ 400-500 eV) detected by ASCA in 1994 
(Turner et al. 1997; Sambruna et al. 1999) suggested 
a standard Seyfert-like scenario with accretion-dominated 
X-rays. BeppoSAX observed NGC 6251 in July 2001 dur- 
ing a high-flux state (the 2 10 keV flux, Fx = 4.7 x 10"^^ 
erg cm~^ s~^, was ~3 times larger than the value mea- 
sured by ASCA 7 years earlier). Based on the absence 
of a prominent Fe Ka line (EW < 100 eV) or other re- 
processing features, Guainazzi et al. (2003) proposed a 
scenario with two main spectral components: a blazar-like 
spectrum dominating the high-flux state and a Seyfert- 
like spectrum emerging during the low-flux state. More 
recently, in March 2002, NGC 6251 was observed with 
XMM-Newton when the 2-10 keV flux was ~15% lower 
than the BeppoSAX value. Despite the relatively high 
flux, the XMM-Newton spectral results seem to support 
the picture emerged from the ASCA observation with the 
presence of a prominent (EW~220 eV) and possibly broad 
{a > 0.3 keV) Fe Ka line (Gliozzi et al. 2004). However, 
the existence of a broad Fe Ka line in the XMM-Newton 
spectrum is still a matter of debate (see Evans et al. 2005 
and Gonzalcz-Martfn et al. 2006 for a discording and a 
supporting view, respectively). 

The controversial results derived from previous X-ray 
studies highlight the impossibility of firmly determining 
the origin of the X-rays in NGC 6251 based solely on 
time-averaged spectral results. Here, we present the re- 
sults from a systematic study of the long-term X-ray flux 
and spectral variability of NGC 6251 using one year long 
proprietary Rossi X-ray Timing Explorer (RXTE) obser- 
vations in the 2-12 keV range. We use model-independent 
methods and time-resolved spectroscopy to study the X- 
ray temporal and spectral properties of this source. The 
main purpose of this analysis is to shed light on the origin 
of the X-rays and in particular on the role played by a 
jet in the X-rays. Once (if) the jet contribution is prop- 
erly assessed, the physical parameters characterizing the 
accretion process onto the supermassive black hole can be 
better constrained, and hence it is possible to discriminate 
between competing theoretical models for the accretion 
process. 

The outline of the paper is as follows. In § 2 we describe 
the observations and data reduction. The main character- 
istics of the X-ray light curve are described in § 3. In § 4 
we study the X-ray spectral variability of NGC 6251 with 
two model-independent methods. In § 5 we describe the 
results of a time-resolved and flux-selected spectral anal- 
yses. In § 6 we discuss the results and their implications. 
In § 7 we summarize the main results and conclusions. 



2. OBSERVATIONS AND DATA REDUCTION 

We use proprietary RXTE data of NGC 6251 that was 
regularly observed for ^ 1000 s once every 4 days be- 
tween 2005 March 4 and 2006 February 27. The obser- 
vations were carried out with the Proportional Counter 
Array (PGA; Jahoda et al. 1996), and the High-Energy X- 
Ray Timing Experiment (HEXTE; Rotschild et al. 1998) 
on RXTE. Here we will consider only PCA data, because 
the signal-to-noise of the HEXTE data is too low for a 
meaningful analysis. 

The PCA data were screened according to the follow- 
ing acceptance criteria: the satellite was out of the South 
Atlantic Anomaly (SAA) for at least 30 minutes, the 
Earth elevation angle was > 10°, the offset from the 
nominal optical position was < 0°02, and the parameter 
ELECTRO N-2 was < 0.1. The last criterion removes data 
with high particle background rates in the Proportional 
Counter Units (PCUs). The PCA background spectra and 
light curves were determined using the L7 — 240 model 
developed at the RXTE Guest Observer Facihty (GOF) 
and implemented by the program pcabackest v. 3.0. This 
model is appropriate for "faint" sources, i.e., those pro- 
ducing count rates less than 40 s~^ PCU~^. 

All the above tasks were carried out using the FTOOLS 
v.6.2 software package and with the help of the REX script 
provided by the RXTE GOF. Data were initially ex- 
tracted with 16 s time resolution and subsequently re- 
binned at different bin widths depending on the applica- 
tion. The current temporal analysis is restricted to PCA, 
STANDARD-2 mode, 2-12.5 keV, Layer 1 data, because 
that is where the PCA is best calibrated and most sensi- 
tive. PCUs and 2 were turned on throughout the mon- 
itoring campaign. However, since the propane layer on 
PCUO was damaged in May 2000, causing a systematic in- 
crease of the background, we conservatively use only PCU2 
for our analysis (see below). All quoted count rates are 
therefore for one PCU. 

The spectral analysis of PCA data was performed us- 
ing the XSPEC V. 12.3. 1 software package (Arnaud 1996). 
We used PCA response matrices and effective area curves 
created specifically for the individual observations by the 
program pcarsp, taking into account the evolution of the 
detector properties. All the spectra were re-binned so that 
each bin contained enough counts for the statistic to be 
valid. Fits were performed in the energy range 2.5-12.5 
keV, where the signal-to-noise ratio is the highest. 

3. THE X-RAY LIGHT CURVES 

Although NGC 6251 is generally considered an X-ray 
bright source - its average flux is of the order of 4 x 10~^^ 
erg cm~^ s~^ with a corresponding luminosity L2_io keV ^ 
5 X 10*^ erg that is nearly a factor 4 larger than the typ- 
ical values observed in low-power radio galaxies (Donato 
et al. 2004) - it is rather weak for the RXTE capabilities. 
Therefore, before starting a detailed analysis of temporal 
properties, it is necessary to demonstrate that the vari- 
ability observed cannot be ascribed to uncertainties in the 
RXTE background or to other artifacts. To this end, we 
have performed the following test: We have compared the 
background-subtracted light curves obtained using PCU2 
layer 1 and PCU2 layer 3. Since the genuine signal in 
layer 3 is quite small, its light curve can be used as a 
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proxy to check how well the background model works. If 
the latter light curve is significantly variable with a pat- 
tern similar to the one produced using layer I, then the 
variability is simply due to un-modeled variations of the 
background. Conversely, if the PCU2 layer 3 light curve 
does not show any pronounced variability or if the flux 
changes are uncorrelated with those observed in the layer 
1 light curve, we can safely conclude that the variability 
detected in NGC 6251 is real. 

The two light curves in the 2-10 keV range (where the 
background PCA model is better parameterized; see Ja- 
hoda et al. 2006 for more details) are shown in Figure [1] 
A visual inspection of this figure suggests that the vari- 
ability in the layer 1 light curve is much more pronounced 
than the one observed in layer 3. Indeed, the mean count 
rate of layer 3 is consistent with zero, indicating that on 
long timescales the background model works adequately. 
However, statistically speaking both light curves are con- 
sidered variable: = 229.5 (layer 1) and 128.3 (layer 
3) for 78 degrees of freedom (hereafter dof). On the 
other hand, a formal analysis based on the excess variance 
(i.e. the variance corrected for statistical errors), indicates 
that variability associated with the layer 1 light curve is 
more than one order of magnitude larger than that asso- 
ciated with layer 3 {a^s = 4 x 10~^ and 3 x 10"'' s~^, 
respectively). Further support to the fact that the vari- 
ability associated with layer 1 is genuine comes from a 
correlation analysis: When layer 3 is plotted versus the 
layer 1 count rate (see Figure [2|), no correlation is ob- 
served, as confirmed by a least square linear fit analysis, 
y = 0.002(±0.015) - 0.005(±0.048)a;. We therefore con- 
clude that most of the count rate changes observed in layer 
1 are associated with genuine intrinsic variations of the X- 
ray source in NGC 6251. 

3.1. PCUO versus PCU2 

In order to maximize the signal-to- noise (S/N) ratio of 
the light curves, one can combine the 2 PCUs at work 
during the monitoring campaign (i.e., PCUO and PCU2), 
provided that they are consistent with each other. In par- 
ticular, it is necessary to test whether PCUO, partially 
damaged since May 2000, is compatible with PCU2. Fig- 
ure [3] shows the background-subtracted, 2-12 keV light 
curves of PCU2 (top panel), PCUO (middle panel), each 
with superimposed the light curve of the other PCU (rep- 
resented by continuous lines in top and middle panels), 
and their ratio (bottom panel) with a time bin of 4 days. 
A visual inspection of this figure suggests that the PCUO 
and PCU2 light curves are not fully consistent with each 
other. Indeed, a formal check based on a test indicates 
that the ratio PCU0/PCU2 is not consistent with the hy- 
pothesis of constancy (x^ = 160.5 for 54 dof). Since the 
propane layer on PCUO was damaged a few years ago, we 
decided to work with the PCU2 data only. For complete- 
ness, we have also performed the same data analysis on 
the PCUO data, and the results are in general agreement 
with those from the PCU2. 

Whereas for the temporal analysis the PCUO data need 
to be necessarily excluded because they appear to be 
inconsistent with the PCU2 data, PCUO spectral data 
should not be a priori ruled out for an analysis averaged 
over two-month periods. Unfortunately, also the 2-month 



averaged PCUO spectra appear to be at odds with the 
PCU2 spectral results, with unphysically large changes for 
the photon index (with a seesaw trend between F ~ 2 and 
F ~ 3) and very large uncertainties. Therefore, also for 
the spectral analysis (see §5) we use only PCU2 data. 

Using data from PCU2 only, we constructed background 
subtracted light curves in two energy bands, namely the 
"soft" (2.5-5 keV) and the "hard" (5-12 keV) band. We 
show them in Figure 21 together with the hardness ratio 
(hard/soft) plotted versus time. Time bins are 4 days. 
The soft and hard band light curves are significantly vari- 
able (x^ = 113.1/78 and x^ = 175.6/78, respectively) and 
appear to be qualitatively similar. This is formally con- 
firmed by the hardness ratio light curve, which appears to 
be constant (x^ = 70.7/78). 

We conclude that the flux variations observed in 
NGC 6251 on timescales of weeks-months are not asso- 
ciated with spectral variations. 

4. SPECTRAL VARIABILITY: MODEL-INDEPENDENT 
ANALYSIS 

For this study we use simple methods such as hardness 
ratio versus count rate plots and the fractional variability 
versus energy plots. These can provide useful information 
without any a priori assumption regarding the shape of 
the X-ray continuum spectrum. Thus, the results from 
the study of these plots can be considered as "model- 
independent" . 

Figure [5] shows the Hard/Soft X-ray color (5-12 
keV)/(2.5-5 keV) plotted versus the total count rate (2.5- 
12 keV) for un-binned (smaller, gray symbols) and binned 
data (larger, darker symbols), respectively. A visual in- 
spection of this figure suggests the presence of a positive 
trend between HR and the count rate (i.e., the spectrum 
hardens as the flux increases), although within a large scat- 
ter. This positive trend is apparently confirmed by the 
linear fit of the data, y = (0.2 ± 0.2) + (2.3 ± 0.5)2: with 

/dof — 43.9/77, which suggests the presence of a posi- 
tive correlation at 4.6(t confidence level. This analysis was 
performed using the routine fitexy (Press et al. 1997) 
that accounts for the errors not only on the y-axis but 
along the x-axis as well. 

Better insight into the presence of correlations between 
HR and total count rate can be obtained by investigat- 
ing the binned data. To this end, we have binned the 
data, shown in Fig. [5l using count rate bins of fixed size 
(0.05 counts/s), and computed the weighted mean (along 
both the y and x-axis) of all the points which fall into a 
bin. The error on the weighted mean is computed follow- 
ing Bevington's prescriptions (Bevington 1969). Only bins 
with at least 5 data points have been plotted in Figure O 
and considered for the linear fit. In this case, the best lin- 
ear fit, y= (0.4±0.2)-h(1.6±0.6) with x^/dof = 1.44/4, 
(shown in Fig. [5] with its la uncertainties), indicates that, 
if a few outliers with large count rate (and error-bar) are 
neglected, the significance of a positive correlation reduces 
to 2.67f7. This is slightly lower than the 3cr level, which 
is generally accepted for significant variations, but it still 
implies a positive correlation at a confidence level of 
99%. 

Another simple way to quantify the variability proper- 
ties of NGC 6251 without considering the time ordering 
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of the values in the light curves, is based on the fractional 
variability parameter i^var ■ This is a commonly used mea- 
sure of the intrinsic variability amplitude relative to the 
mean count rate, corrected for the effect of random errors, 

-Tvar - 

where is the variance, (r) the unweighted mean count 
rate, and the mean square value of the uncertainties 
associated with each individual count rate. The error on 
Fvar has been estimated following Vaughan et al. (2003). 
We computed Fvai- on selected energy bands, with mean 
count rates similar and sufficiently high. 

The results are plotted in Figure[S]and suggest that, over 
the 3-8 keV energy range, a weak positive trend seems to 
be present. However, the large uncertainties on i^var due 
to the low count rate in the narrow energy-selected bands 
hamper this kind of analysis. Indeed, statistically speak- 
ing, the trend shown in Figure [5] is consistent with the hy- 
pothesis of Fvar being constant in the energy range probed 
by RXTE. If only two broader energy bands are used for 
this analysis, the resulting fractional variability in the hard 
energy band 5-12 keV, -Fvar,hard = (20±3)% appears to be 
larger than the value, Fvar,soft = (10±5)%, obtained in the 
soft band 2.5-5 keV. However, the i^var,hard — ^var.soft = 
(10 ± 6)% difference is significant at the 1.1a level only. 

In summary, the results from the study of the HR — ct 
plots show evidence for a positive correlation between 
hardness ratio and total count rate. Similarly, an anal- 
ysis of the fractional variability suggests that Fvar is more 
pronounced in the hard energy band than in the soft one. 
However, due to the rather limited signal-to-noise ratio of 
our data these results are significant at just the 2.7 and 
1.7a level, respectively. 

5. SPECTRAL ANALYSIS 

5.1. Time-resolved Spectroscopy 

Given that the data consist of short snapshots span- 
ning a long temporal baseline, they are in principle well 
suited for monitoring the spectral variability of NGC 6251. 
However, due to the limited S/N of the data, the spectral 
slope r, measured from spectra of individual observations, 
cannot be adequately constrained and hence the spectral 
variability cannot be investigated in such a way. Indeed, 
if we plot the values of T versus time, no variations are 
detected due to the large errors on F. This is fully con- 
sistent with the results from the iJi?-time plot described 
in §4. In order to increase the S/N and investigate the 
presence of possible spectral variations, we use spectra av- 
eraged over two- month intervals. This choice is a trade-off 
between the necessity of accumulate sufficient counts for a 
reliable spectral analysis and the need to use limited tem- 
poral intervals to minimize the effects of the slow drift in 
the detector gain. 

We fitted each two-month spectrum with a simple 
power-law (PL) model absorbed by Galactic Nn (5.65 x 
10^° cm~^). The model fits all the data reasonably well, 
as indicated by Fig. [7] that shows a typical spectrum fit- 
ted with a simple power law. The best-fit results are 
listed in Table 1 and can be summarized as follows. A 
simple PL model provides an acceptable parametrization 



for all spectra. The photon indices are all rather steep 
(F ~ 2.5) and consistent with each other within the er- 
rors. In other words, our results suggest that the source's 
spectrum does not vary significantly on timescales longer 
than two months. The weighted spectral slope mean is 
2.5 ± 0.1. Adding a Gaussian line at 6.4 keV to the PL 
continuum model does not improve the fit significantly in 
any of the six spectra, but the 90% confidence upper limits 
on the equivalent width are relatively high (EW~200-700 
eV). 

In order to better constrain F and the line EW, we have 
tried a comparison between the first and the last 6 months 
of monitoring campaign, by fitting together 3 two-month 
spectra at a time. The results - F2-10 koV = 4.3 x 10"-^^ 
erg cm-2 s-\ F = 2.53 ± 0.18, EW < 212 eV during 
the first half, and F2-10 kcV = 4.5 x 10"^^ erg cm~^ s~^, 
F = 2.46 ± 0.17, EW < 144 eV during the second half - 
are fully consistent with each other. 

Since there are no indications for long-term spectral 
variability, we have fitted the 6 two-month spectra to- 
gether. This yielded F = 2.5 ± 0.1 and a significantly 
smaller upper limit on the line strength: EW<104 eV. In 
the 2-10 keV energy band, we obtained an average ab- 
sorbed fiux of 4.4x 10""'^^ erg cm~^ s""'^, and a correspond- 
ing intrinsic luminosity of 6.3 x 10*^ erg s~^, assuming 
Hq = 71 km s~i Mpc-i, f^A = 0.73 and = 0.27 (Ben- 
net et al. 2003). 

5.2. Flux- selected Spectroscopy 

In order to verify the presence of a direct correlation 
between HR and count rate derived in §4, we performed a 
flux-selected spectral analysis. To this end we divided the 
94 individual spectra into 5 bins according to their average 
count rate (namely, < 0.25, ,0.25-0.30, 0.30-0.35, 0.35- 
0.40, and > 0.40 c/s). In each bin, the individual spec- 
tra were fitted simultaneously with a PL model (absorbed 
by Galactic iVn) with their photon indices linked together 
and their respective normalizations free to vary. 

The results are summarized in Table 2, where the re- 
ported errors on F and flux are respectively la and a/ VN, 
with N being the number of individual spectra per bin. 
The upper limits on EW correspond to the 90% confi- 
dence level. Table 2 indicates that the 2.5-12.5 keV spec- 
tra harden as the average fiux increases. This is clearly 
shown in Figure [51 where the values of F for each count 
rate bin have been plotted against their respective flux 
values. The dashed line represents the best linear fit, 
y = 3.6 ± 0.3 - (0.24 ± 0.08)a;, which reveals that the 
inverse correlation is significant at 3-cr level. 

For comparison, in Fig. [5] we have also plotted the val- 
ues corresponding to the ASCA, BeppoSAX, and XMM- 
Newton observations. To this end, we have used PIMMS to 
convert the observed flux into the RXTE energy range, 
assuming the best-fit spectral parameters reported in 
Guainazzi et al. (2003) for BeppoSAX and ASCA, and 
Gliozzi et al. (2004) for XMM-Newton. The ASCA, Bep- 
poSAX, and XMM-Newton photon indices seem to follow 
the same inverse trend shown by RXTE, becoming flat- 
ter as the flux increases. However, their values appear 
to be significantly smaller than those obtained from the 
flux-selected RXTE spectral analysis. This apparent dis- 
crepancy may be probably reconciled (at least for the Bep- 
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poSAX and XMM-Newton data) by bearing in mind that: 
1) Past studies have shown that photon indices measured 
with the RXTE PC A are systematically steeper than those 
measured by other X-ray satellites (e.g., Yaqoob 2003). 2) 
The error-bars for the RXTE data in Fig. [5] are at Icr 
level. On the other hand, it appears more difficult to rec- 
oncile the ASCA value with the extrapolation of the RXTE 
best-fit trend, suggesting that the source was in a different 
physical state (for example, it lacked a jet contribution). 

We have also added a Gaussian line to the PL continuum 
(see Table 2 last column) . Only for the bin with the lowest 
flux, the line is very marginally significant (the value of 
decreases by 1.2 for one additional dof) with EW~400 eV 
and large uncertainties. For bins with larger flux values, 
only upper limits on EW can be derived. These findings 
indicate that a Fe Ka line is never statistically required, 
although its presence cannot be completely ruled out when 
the source flux is at the lowest level measured by RXTE. 

6. DISCUSSION 

We have undertaken the first X-ray monitoring study 
of the FR I galaxy NGC 6251, investigating the tempo- 
ral and spectral variability as well as time-averaged and 
flux-selected spectral results. 

By comparing the X-ray fluxes measured by differ- 
ent satellites over nearly a decade of observations, it is 
clear that NGC 6251 shows large flux changes on long 
timescales. For example, in October 1994 ASCA mea- 
sured a flux of ~ 1.4 X 10~^^ erg cm^^ s^^ in the 2-10 
keV energy band, whereas BeppoSAX observed a flux of 
- 4.7 X 10-12 gj.g cm-2 s-i in July 2001, and XMM- 
Newton of ~ 4 X 10-12 erg cm-^ s-^ in March 2002. The 
high throughput of the EPIC pn camera aboard XMM- 
Newton also revealed the presence of low-amplitude flux 
changes (of the order of ~5%) on timescales of a few ks. 

The variability behavior detected in the current RXTE 
observations is in agreement with the sparse evidence that 
had been gathered in the previous years, and show con- 
clusively that NGC 6251 is a persistently variable source 
in X-rays. The frequent RXTE observations, spread over 
a period of one year, indicate that this radio galaxy is 
characterized by persistent variability in the total (2.5-12 
keV), soft (2.5-5 keV), and hard (5-12 keV) energy bands. 
Throughout the observation, the flux appears to randomly 
change by a factor of ^2 on timescales of a few weeks. 

Previous X-ray studies, mostly based on the time- 
averaged spectral analysis, have revealed a dual behavior 
for NGC 6251: sometimes the source appears to behave 
like a Seyfert galaxy with a strong Fe Ka line (e.g.. Turner 
et al. 1997; Sambruna et al. 1999; Gliozzi et al. 2004), 
whereas in other occasions NGC 6251 seems to be more 
consistent with a blazar-like behavior, showing a feature- 
less X-ray spectrum (e.g., Guainazzi et al. 2003). The 
densely sampled, year-long RXTE observations and the 
investigation of the flux and spectral variability properties 
offer an alternative and model-independent way to shed 
some light on this issue. 

6.1. Evidence from the flux variability properties 

Unlike the brightest blazars frequently monitored by 
RXTE (e.g., Kataoka et al. 2001, Cui 2004, Xue & Cui 
2005), NGC 6251 does not show any prominent flare on 



any observable timescale, nor does it show the large am- 
plitude variability typically observed in several Seyfert- 
like objects during yearly-long monitoring campaigns (e.g., 
Markowitz & Edelson 2001). 

On one hand, the apparent inconsistency with the large 
variability observed in Seyfert galaxies can be explained 
by the lower values of black hole mass in the latter ob- 
jects, which is typically one order of magnitude lower than 
NGC 6251. On the other hand, the lack of prominent flares 
in the NGC6251 light curve, which instead characterize the 
blazar light curves, can be understood by keeping in mind 
that the blazars monitored by RXTE are the brightest 
members of this AGN class and that the observations are 
often triggered only during their flaring activity. Nonethe- 
less, blazar monitoring campaigns with baselines covering 
several years reveal that also the brightest blazars alter- 
nate prominent flaring activity with "quiescent periods" 
that are characterized by moderate flux variations. In- 
deed, prolonged periods of moderate variability have been 
detected in several blazars (e.g., Bottcher et al. 2005; 
Marscher 2006). 

A model-independent study of the RXTE variability 
properties of the prototypical blazar Mrk 501, similar to 
the one performed on NGC 6251 in this work, offers the 
best opportunity for a more quantitative comparison of 
NGC 6251 with a typical blazar behavior (Gliozzi et al. 
2006). For instance, Mrk 501, which showed a large out- 
burst in 1997, underwent a progressive decrease of its ac- 
tivity in the following years, resulting in a lower mean 
count rate accompanied by lower variability. Speciflcally, 
in 1999, when Mrk 501 reached a minimal flux value, the 
light curve was characterized by variations of the order of 
~2 on timescales of 20-30 days, which are fully consistent 
with those detected in NGC 6251. 

Before proceeding further, it must be pointed out that 
the blazars regularly monitored over the past years with 
RXTE (including Mrk 501) are basically all High-peaked 
blazars (HBLs), with the synchrotron component peak- 
ing in the X-ray range and a second spectral component, 
generally attributed to inverse Compton scattering, peak- 
ing at TeV energies. On the other hand, in the blazar 
framework, the broadband SED of NGC 6251 is consistent 
Low-peaked blazars (LBLs), with the inverse Compton 
component peaking in the X-rays (Chiaberge et al. 2003; 
Guainazzi et al. 2003). As a consequence, a formal com- 
parison of the X-ray properties of NGC 6251 should be in 
principle performed using the TeV properties of Mrk 501. 
However, detailed studies of X-ray and TeV emissions in 
HBLs have demonstrated the existence of a tight correla- 
tion between these energy bands, indicating that the X-ray 
and TeV radiation follow the same variability trend (e.g., 
Fossati et al. 2004; Ghozzi et al. 2006). Further sup- 
port to this conclusion comes from recent investigations of 
the TeV properties of Mrk 501 and Mrk 421 carried out 
with the MAGIC and Whipple telescopes (Blazejowki et 
al. 2006; Albert et al. 2007). These studies demonstrate 
that also in the 7-ray energy band HBLs show the typi- 
cal blazar spectral variability behavior observed in X-rays, 
that is a spectral hardening when the source brightens and 
a fractional variability more pronounced at higher energies 
(e.g., Zhang et al. 1999; Fossati et al. 2000; Gliozzi et al. 
2006; Rebillot et al. 2006 and references therein). We 
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therefore conclude that it is appropriate to compare the 
X-ray properties of Mrk 501 with those of NGC 6251. 

Another relevant finding from the flux variability anal- 
ysis of NGC 6251 is that the fractional variability appears 
to be more pronounced in the hard than in the soft en- 
ergy band: AFvar = -FVar.hard - ^Var.soft = (10 ±6)%. This 
result can be quantitatively compared with similar stud- 
ies carried out on Scyfcrt galaxies monitored with RXTE 
for several months or years (Markowitz & Edelson 2004). 
Speciflcally, all the Seyfert galaxies with long i?^?^? mon- 
itoring yield negative values of AFvar ranging between -3.5 
and -24.8, with a mean value of — 10.3± 1.5. On the other 
hand, a similar study, carried out using RXTE monitoring 
data of the blazar Mrk 501 between 1997 and 2000, yielded 
positive values of AFyar ranging between 5 and 42, with a 
mean value of 15 ± 1. 

We can therefore conclude that the limited amplitude 
of the NGC 6251 flux variability can be explained equally 
well in both Seyfert-like and blazar-like scenarios. How- 
ever, the fractional variability behavior is inconsistent with 
the trend generally observed in Scyfcrt galaxies, but fully 
consistent with the typical blazar behavior observed in 
Mrk 501. 

6.2. Evidence from the spectral properties 

A comparative analysis of the spectral variability of 
NGC 6251 and the typical behaviors observed in both 
radio-quiet AGNs (where the X-rays are thought to be 
produced by Comptonization in the corona that is closely 
connected with the accretion disk) and radio-loud jet- 
dominated AGNs (whose radiation over the entire energy 
range is ascribed to jet emission) can help us understand 
whether the X-ray radiation from NGC 6251 is dominated 
by jet or accretion-related emission. 

A flux-selectcd spectral analysis (in agreement with the 
HR — ct plot) has shown that the X-ray spectrum of 
NGC 6251 hardens as the source brightens, following the 
linear correlation F oc — 0.24(±0.08) x Fx- A direct com- 
parison of this result with a similar spectral study, carried 
out on 4 Seyfert galaxies monitored with RXTE, indicate 
that the latter have always positive slopes in the F — Fx 
plane ranging between 0.05 and 0.15 with a mean value 
of 0.09 ± 0.01 (Papadakis et al. 2002). As a consequence, 
the spectral variability behavior of NGC 6251 appears to 
be inconsistent with the typical Scyfcrt-likc trend. On the 
other hand, the F — Fx slope measured for Mrk 501 during 
the weakly variable period of 2000, —0.37 ± 0.06, appears 
to be fully consistent with the behavior of NGC 6251 dur- 
ing the RXTE monitoring campaign. 

The existence of a Seyfert-like component suggested by 
previous X-ray studies was essentially based on the pres- 
ence of a strong Fe Ka line. Unfortunately, the low S/N 
spectra obtained during the RXTE monitoring coupled 
with the poor spectral resolution of the RXTE PCA ham- 
pers a detailed investigation of this issue. Indeed, if the 
2-month averaged PCA spectra are fitted with a model 
including a power law and a Gaussian line with spectral 
parameters fixed at the best-fit values obtained during the 
ASCA observation (we conservatively assumed the best-fit 
parameters reported in Guainazzi ct al. 2003), the results 
are statistically indistinguishable from the fits obtained 
using a simple power law. This indicates that RXTE is 



unable to confirm or refute the presence of a Fe Ka line, 
and suggests that only a relatively long exposure of XMM- 
Newton with its superior capabilities is able to detect the 
presence of the line when the average fiux of NGC 6251 is 
relatively high. 

6.3. The origin of X-rays in NGC 6251 

The primary goal of this work is to investigate the ori- 
gin of the X-rays in NGC 6251 and in particular to assess 
the role played by the putative jet. At first sight, the 
possibility that the jet may dominate the radiative out- 
put of a radio galaxy may be surprising, and even more 
so for NGC6251, which has a Mpc radio jet (and hence 
appears forming a large viewing angle) and Seyfert-like 
emission during the low flux state. However, under spe- 
cific circumstances, the jet-dominance hypothesis becomes 
plausible. This is the case when the base of the jet is not 
well collimated and the X-rays are produced by a part of 
the outfiowing material that points towards the observer. 
Indeed, this is the framework proposed to explain the TeV 
emission detected in another radio galaxy M87 (Aharonian 
et al. 2006). Alternatively, if the base of the jet is tilted 
with respect to the large scale jet and forms a small view- 
ing angle, the possible dominance of the jet in the 2-12 
keV range can be naturally explained. In fact, this is the 
scenario put forward by Jones & Wehrle (2002) to explain 
the large jct/countcrjet brightness ratio on parsec scales 
inferred for NGC 6251 from VLBA observations. 

By comparing the spectral variability properties of 
NGC 6251 with those of radio-quiet AGNs and blazar 
objects, we find that they are certainly inconsistent with 
a typical Seyfert-like behavior but fully consistent with 
blazars. Combining pieces of information from the model- 
independent analysis with the findings from the flux- 
selected spectral analysis, we are led to the conclusion 
that, during the RXTE monitoring campaign, the bulk 
of the hard X-ray radiation from NGC 6251 was domi- 
nated by the emission from the unresolved base of the jet. 
Nonetheless, the presence of a disk-corona component, de- 
tected in previous observations with ASCA and XMM- 
Newton, cannot be ruled out; the upper limits measured 
on the equivalent width of the Fe Ka line are indeed fully 
consistent with the values measured by ASCA and XMM- 
Newton, but the low S /N of the spectra hampers a more 
quantitative analysis. 

In conclusion, we can try to exploit the main results 
from this work (i.e., the jet dominance in the high flux 
state, with a possible contribution from a Seyfert-like com- 
ponent emerging at low flux values) to derive some con- 
straints on the accretion process at work in NGC 6251. 
Assuming that an accretion-related component is always 
present in NGC 6251 and dominates during the low flux 
state, we can use the average flux measured in the low 
count rate bin (sec Table 2) to compare the correspond- 
ing 2-10 keV Imninosity - Lx = 6.6 x 10^^ erg s~^- 
to the Eddington value readily derived from the black 
hole mass estimate. The relatively high value derived, 
Lx/LEdd > 5 X 10"'"', confirms that NGC 6251 is a bright 
FR I galaxy that is close to the FR I/FR II dividing line 
in terms of power of the central engine. However, it is not 
possible to derive any firm conclusion on the nature of the 
accretion process given the unknown contribution of the 
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jet in the low-fliix state. 

7. SUMMARY AND CONCLUSIONS 

We have used data from a year-long RXTE monitoring 
campaign to study the spectral variability of NGC 6251 

following model-independent and spectral model fitting 
methods. The main results can be summarized as follows: 

• Throughout the monitoring campaign and espe- 
cially during the last 4 months, NGC 6251 was in 
a relatively high-flux state, with values of the 2-10 
keV absorbed flux comparable to that observed by 
BeppoSAX in 2001. 

• The light curves show persistent variations by a fac- 
tor of 2 on timescales weeks/months in the total 
(2.5-12 keV), soft (2.5-5 keV), and hard (5-12 keV) 
energy bands. 

• The fractional variability, computed over the soft 
and hard energy bands, reveals an enhanced vari- 
ability at harder energies (Fvar ^20%) compared to 
the soft band (F^ar --10%). 

• There is evidence of a positive trend in the HR — ct 
plot (or, analogously, of a negative trend in the F- 



flux plot); in other words the spectrum hardens as 
the flux increases. 

• The 2-month averaged spectra are well fitted by a 
power-law model, with (F) ~ 2.5 and no indication 
for an Fe Ka line. Combining all the 2-month av- 
eraged spectra yielded EW< 104 eV. Only for the 
lowest flux-sclcctcd spectrum, there is marginal ev- 
idence for a Fc Ka line. However, the low S/N does 
not allows one to put reliable constraints on the 
putative Seyfert-like component. 

The inconsistency of the spectral variability behavior 

of NGC 6251 with the typical trend observed in Seyfert 
galaxies and the similarity with blazars lead support to a 
jet-dominated scenario. However, based on the RXTE ob- 
servations, a substantial contribution from a disk-corona 
system cannot be ruled out. 

We thank the referee for the comments and suggestions 
that improved the clarity of the paper. MG acknowledges 
support by the RXTE Guest Investigator Program under 
NASA grant 200858. Funds from the NASA LTSA grant 
NAG5-10708 are also gratefully acknowledged. 
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Fig. 1. — Top panel: RXTE PCA light curve in the 2-10 keV energy band, using PCU2 layer 1. Bottom panel: PCU2 layer 3 light curve. 
Time bins are 4 days. The dashed lines are the average count rate level. 




Fig. 2. — Layer 3 plotted versus the layer 1 count rate. The dashed line represents the best linear fit. 
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Fig. 3. — Top panel: PCU2 light curve of NGC 6251 in the 2-10 keV range; the solid line represents the PCUO light curve. Middle panel: 
PCUO light curve of NGC 6251 in the 2-12 keV range; the solid line represents the PCU2 light curve. Bottom panel: PCU0/PCU2 light 
curve; the dashed line represents the average value of the ratio PCU0/PCU2. Time bins are 4 days. 
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Fig. 4.— Top panel: Soft (2.5-5 keV) light curve of NGC 6251 using PCU2 data only. Middle panel: Hard (5-12 keV) light curve. Bottom 
panel: Hardness Hard/Soft ratio light curve; the solid line represents the average value of the hardness ratio. Time bins are 4 days. 
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Fig. 5. — Hardness ratio (5-15 keV/2.5-5 keV) plotted versus the total count rate. The gray (light blue in color) small symbols correspond 
to the un-binned data points, whereas the darker, larger symbols refer to the weighted mean of data points in bins of fixed size of 0.05 c/s 
along the x-axis. The continuous line represents the best linear fit from the binned data, whereas the dashed lines are the ztlcr departure 
from the best fit. 
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Fig. 6. — Fractional variability amplitude as a function of energy for NGC 6251. The error-bars along the x axis simply represent the 
energy band width. The error bars along the y axis are computed following Vaughan et al. 2003. 
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Fig. 7. — PCA spectrum of NGC 6251 during the third 2-month interval obtained using PCU2 data only. The model used is a simple power 
law absorbed by Galactic Nn. 
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Fig. 8. — r, obtained from the spectral analysis of flux-selected intervals, plotted against the average flux in the 2.5-12.5 keV band. The 
error-bars on F are l-tr errors. The dashed line represents the best linear fit y = 3.6 ± 0.3 — (0.24 ± 0.08)a;. For comparison, we have also 
plotted the values corresponding to the ASCA, BeppoSAX, and XMM-Newton observations. To this end, we have used PIMMS to convert the 
observed flux into the RXTE energy range, assuming the best-fit spectral parameters reported in Guainazzi et al. (2003) and Gliozzi et al. 
(2004). 
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Table 1 

Time-resolved Spectral results 



Interval 


Xred 


^2-10 keV 


r 


EWpeKa 

(eV) 


1 


0.60 


4.1 


2 6+°-'' 

■^■'^-0.4 


< 366 


2 


0.53 


4.5 


'^•"-0.3 


< 466 


3 


0.45 


4.4 


9 t;+U.3 
''•^-0.4 


< 444 


4 


0.81 


4.0 


2 4+^'-'* 

^■^-0.4 


< 756 


5 


1.00 


5.2 


9 7+0.3 
^■'-0.3 


< 181 


6 


0.59 


4.8 


9 4+0.3 


< 246 



Note: All errors and upper limits refer to 90% confidence levels. 



Table 2 

Flux-selected Spectral results 



Individual Spectra 


Count Rate 


Xred 


^2. 5-12. 5 keV 


r 


EWpcKa 




(s-i) 




(10~^^erg cm~^ s~^) 




(eV) 


31 


< 0.25 


0.53 


2.5 ±0.1 


2.9 ±0.2 


< 903 


26 


0.25 - 0.30 


0.55 


3.2 ±0.1 


2.7 ±0.2 


< 540 


16 


0.30-0.35 


0.50 


3.8 ±0.1 


2.5 ±0.2 


< 205 


12 


0.35 - 0.40 


0.59 


4.4 ±0.1 


2.8 ±0.2 


< 378 


9 


> 0.40 


0.60 


5.3 ±0.3 


2.2 ±0.2 


< 305 



Note: The errors on T are la (68% confidence level), whereas the upper limits on EW are 90% confidence levels. 



